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Abstract 

The numerical solution of the Navier-Stokes equations 
about three-dimensional configurations is reviewed. For- 
mulational and computational requirements for the various 
Navier-Stokes approaches are examined for typical prob- 
lems including the viscous flow field solution about a com- 
plete aerospace vehicle. Recent computed results, with 
experimental comparisons when available, are presented 
to highlight the presentation. The future of Navier-Stokes 
applications in three-dimensions is seen to be rapidly ex- 
panding across a broad front including internal and exter- 
nal flows, and flows across the entire speed regime from 
incompressible to hypersonic applications. Prospects for 
the future are described and recommendations for areas of 
concentrated research are indicated. 

Introduction 

The fields of Computational Fluid Dynamics (CFD) 
and computer hardware design have been rapidly advanc- 
ing over the last decade. It is now possible to demonstrate 
the importance of CFD in the field of aircraft design and 
performance prediction for general flow conditions includ- 
ing reasonably complete viscous effects. The problem of 
solving the Navier-Stokes equations about complex three- 
dimensional geometries has been undertaken and signifi- 
cant results have already been produced. The purpose of 
this paper is to review the current status of CFD applica- 
tions in this area and to give prospects for the future espe- 
cially for those applications involving complex geometries. 
Computational requirements for producing these results 
both now and in the future will be addressed. 

The first pioneering calculations with the Navier- 
Stokes equations generally involved turbulent solutions 
with simple algebraic turbulence models (or laminar solu- 
tions with no turbulence models) in two dimensions with 
generally simple geometric shapes. These aspects were re- 
quired to eliminate the complexity of turbulence and to 
reduce computational requirements to fit the computa- 
tional facilities that existed in the early and mid 1970s. 
A good review of early work in this area is given by Peyret 
and Viviand. 1 Several years after the first two-dimensional 
Navier-Stokes solutions appeared, three-dimensional sim- 
ulations were presented, again using algebraic turbulence 
models (or a laminar flow assumption) and simple geome- 
tries. Since then improvements in computers and algo- 
rithms have allowed tremendous advances in CFD across 
a broad front of different applications. 

The rate of increase in the number and size of three- 
dimensional Navier-Stokes simulations is displayed in Fig. 
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1. This figure shows grid size for a representative sample of 
three-dimensional Navier-Stokes simulations plotted ver- 
sus the year the results were published. Several observa- 
tions are immediately obvious. The first three-dimensional 
simulations did not start appearing until the mid- 1970’s 
and since then the growth in both number and size of 
simulations has been almost exponential. This trend is 
due to improved algorithms, faster computers with larger 
main memories and to the availability of more computer 
time. The approximate cost of these computations can be 
inferred from Fig. 2 where the number of operations is 
plotted versus grid size. Note that the speed regime asso- 
ciated with each computation ranging from incompressible 
to hypersonic is also indicated. The solid line is an esti- 
mate of projected growth that appeared in Ref. 2 in 1983. 
This estimate was based on only the studies available in 
the pre-1983 time frame, and yet the curve represents a 
reasonablly valid trend. Actually, the number of opera- 
tions required to obtain a solution for the finer grids is 
somewhat underpredicted by this estimate. This is proba- 
bly due to additional “stiffness” associated with finer grids 
and the fact that these solutions are newer and therefore 
probably less efficient. 

This survey is divided^nto several sections listed as 
follows: CFD Formulations, where the major governing 
equations used in CFD are described with special em- 
phasis on the Navier-Stokes formulation; Previous Survey 
Work, where past surveys are listed; Status of Navier- 
Stokes CFD, where Navier-Stokes CFD is surveyed includ- 
ing the areas of incompressible flow, transonic flow, super- 
sonic/hypersonic flow and high-alpha/vortex-dominated 
flow; and Concluding Remarks. Because of the breadth 
associated with this review and the Navier-Stokes CFD 



YEAR 

Fig. 1. Growth in the number and size of three-dimension- 
al Navier-Stokes simulations using the time-asym- 
ptotic approach to obtain steady solutions. 
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Fig. 2. Estimated floating-point operations required for a 
converged solution as a function of grid size for 
selected three-dimensional time-asymptotic RANS 
solutions. 

field in particular, no claim of completeness is made nor 
should one be assumed. The first item for discussion is 
how the Navier-Stokes formulation, and in particular, the 
Reynolds-Averaged Navier-Stokes (RANS) formulation, 
fits into all the other CFD formulations. 


CFD Formulations 

The aforementioned improvements in computer exe- 
cution speed and main memory capacity achieved in the 
past several years has made it possible to dramatically 
expand the field of CFD. This has resulted in significant 
progress up the hierarchal “ladder” of CFD formulation 
complexity as shown in Fig. 3. Panel method solutions, 
which effectively solve the linear Laplace equation using 
the principle of superposition, have been studied thor- 
oughly during the past three decades. As a result, panel 
method computer codes are well established within indus- 
try as efficient and reliable design tools for aircraft in the 
subsonic and supersonic speed regimes. 


The next rung on the CFD formulational ladder is oc- 
cupied by nonlinear potential formulations including the 
transonic small-disturbance and full potential equations. 
Both of these formulations have been extensively studied 
for the past 15 years and except for some grid genera- 
tion limitations associated with complex geometry appli- 
cations, these methods are well developed and heavily uti- 
lized within the aircraft industry design environment. Of 
course, these potential methods extend the range of appli- 
cability of CFD to the transonic speed regime and there- 
fore have the capability of capturing weak shock waves and 
wave-drag levels. Appropriate reviews of this area can be 
found in Refs. 3 and 4. 

The Euler equation formulation, which occupies the 
next rung on the CFD formulational ladder, is currently 
at the center of research in the Aerospace CFD commu- 
nity. New algorithms and applications in three dimensions 
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Fig. 3. Computational Fluid Dynamics formulational lad- 
der. 


have proliferated at a very fast rate in the past live years. 
The various design groups within industry have already 
started using Euler codes for many applications but sig- 
nificant improvements in code and/or computer efficiency 
have to be realized before the Euler formulation will be 
used extensively. The Euler formulation provides better 
accuracy for strong shock calculations and has the abil- 
ity to properly track vorticity gradients. These features 
are important for a variety of advanced flow field situa- 
tions including vortical aerodynamic studies, strong shock 
calculations, high angle-of-attack leading-edge flow sepa- 
ration, high-speed blunt body problems, and blast wave 
studies. 

The highest entry on the CFD formulational ladder 
is occupied by the Navier-Stokes equations. This formula- 
tion is the most complex and represents the smallest level 
of physical approximation. Like the Euler formulation, 
much work has been completed in the three-dimensional 
Navier-Stokes area during the past five years. However, 
the application of typical Navier-Stokes computer codes in 
the industry design environment has been slow to evolve. 
This is largely due to two reasons: First, the Navier-Stokes 
equations are several times more expensive to solve (at 
least) than the Euler equations. This is due to more terms 
in the Navier-Stokes equations; more grid points required 
to resolve a given problem, because in addition to the in- 
viscid region, the boundary layer has to be resolved; and 
in separated flow cases, more iterations are required to 
achieve convergence because of the more complicated vis- 
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cous/inviscid interaction. The second major reason that 
the Navier- Stokes equations have not been used in design 
applications is due to the complications of turbulence. To 
discuss this subject a further breakdown of the Navier- 
Stokes formulation is required. 

The Navier-Stokes equations intimately describe the 
motion of most fluid flow phenomena including large scale 
motions associated with the inviscid features of a flow 
field and small scale motions associated with turbulence. 
The numerical procedure which directly solves the Navier- 
Stokes equations and resolves all flow length scales includ- 
ing the smallest turbulent length scales is called direct sim- 
ulation (DS). When the numerical grid resolution is not 
capable of resolving all length scales problems occur. Be- 
cause for high Reynolds number flow the range of length 
scales is so immense, it is not currently possible to simulate 
this type of flow with schemes based on the DS approach. 
Thus, the first level of approximation to the Navier-Stokes 
equations, called Large Eddy Simulation (LES), is intro- 
duced. In the LES formulation the Navier-Stokes equa- 
tions are solved on a grid which is capable of resolving all 
length scales except the smallest turbulent length scales. 
These length scales are associated with the smallest turbu- 
lent eddies in the flow field. Because these smallest eddies 
tend to be isotropic, they are easier to model than the 
larger eddies which tend to be random and can vary in 
nature for different flow conditions. A good review of the 
DS and LES areas is given in Ref. 5. 

The LES formulation is still quite expensive and is re- 
stricted to problems with limited geometrical complexity 
and to moderate and low Reynolds numbers. This leads to 
the Reynolds-averaged Navier-Stokes (RANS) formulation 
where all turbulent eddies are modeled. The RANS formu- 
lation is derived by decomposing all dependent variables, 
e.g., density, pressure, and velocity components, into two 
parts, a mean part and an unsteady part. The unsteady 
part contains the unsteadiness due to high-frequency tur- 
bulent fluctuations and the mean part contains the aver- 
age local flow value for that variable including any low- 
frequency time variation. After this variable decompo- 
sition process the resulting Navier-Stokes equations are 
time averaged over a time scale which is large in compar- 
ison with turbulent fluctuations but small in comparison 
with mean flow field unsteadiness. As a result of the time- 
averaging process, certain new terms arising from turbu- 
lent eddy correlations remain. These new terms account 
for the mean or average effect of turbulence and must be 
modeled. This allows the numerical solution of the Navier- 
Stokes equations for many applications involving realistic 
Reynolds numbers on grids that will fit on existing super- 
computers. 

One variation of the RANS equations first suggested 
by Baldwin and Lomax 6 and used in a variety of applica- 
tions, is the so-called “thin-layer” formulation. This for- 
mulation requires a no-slip boundary surface to be aligned 
with a constant coordinate surface, which is a typical facet 
in most applications. The thin-layer formulation is ob- 
tained by dropping all viscous terms containing transverse 


spatial derivatives. The justification for this assumption is 
that these terms are typically small in most regions of flow 
and that in regions where they are potentially large, for ex- 
ample, near separation, grid resolution is not generally ad- 
equate to resolve them. In addition, numerical experimen- 
tation using both the full and thin- layer forms has shown 
little or no difference in the computed results, and the thin- 
layer formulation is approximately 30% cheaper. It should 
be stressed that although the thin-layer Navier-Stokes as- 
sumptions are similar to certain classical boundary layer 
assumptions, the resulting equations are quite different. 
All time terms are retained, that is, the thin-layer Navier- 
Stokes equations are still hyperbolic, while the boundary 
layer equations are parabolic. Thus, the thin-layer Navier- 
Stokes formulation has no numerical difficulty with either 
streamwise or cross-flow separations. 

The last Navier-Stokes formulation discussed here is 
the Parabolized Navier-Stokes (PNS) formulation. This 
formulation is derived from the RANS equations and there- 
fore requires the modeling of all turbulent length scales. 
The PNS approach also requires the flow to be supersonic 
in the streamwise direction (although under certain condi- 
tions when the pressure field is specified the freestream can 
be subsonic) and that there be no streamwise separations 
in the boundary layer. In addition, all time-dependent 
terms and streamwise diffusion terms are neglected. Un- 
der these assumptions the Navier-Stokes equations become 
parabolic in the streamwise direction and can be solved 
numerically using a space-marching procedure instead of 
the standard time- asymptotic approach. Thus, the com- 
puter requirements (both memory and cpu time) are much 
less demanding for the PNS approach relative to the time- 
dependent approach. 

A comparison of the estimated memory and execu- 
tion time requirements for several different formulations 
is given in Fig. 4. An additional horizontal scale indi- 
cates the amount of execution time required by a the- 
oretical gigaflop computer, which should be available in 
the near-term future. This figure was constructed from 
information given in Refs. 7-11. The Euler formulation 
and three of the Navier-Stokes formulations just discussed 
(DS, LES, and RANS) are all included. Note that sev- 
eral variations of the RANS formulation (steady ideal gas, 

steady nonequilibrium real gas, and unsteady ideal gas) 
are presented. Including nonequilibrium real gas effects 
associated with hypersonic flight or the unsteadiness ef- 
fects associated with (for example) high-alpha flight, may 
increase the computational costs of a typical steady RANS 
calculation by an order of magnitude. This additional ex- 
pense while prohibitive today may not be prohibitive in 
the near-term future. 

The prospects of computing three-dimensional flows 
at realistic flight Reynolds numbers using the LES formu- 
lation in the near future do not seem bright. As indicated 
in Fig. 4 the cpu time for such a calculation on a theoreti- 
cal gigaflop computer would be on the order of one month. 
The DS picture is even more dismal with cpu times for a 
single three-dimensional calculation requiring on the order 


283 




1 1 1 1 1 1 

1 sec 1 min 1 hr 1 day 1 mon 1 yr 
CPU TIME ON A GIGAFLOP COMPUTER 

Fig. 4. Estimated memory and execution time require- 
ments for several different CFD formulations. 

of years. Hence, it seems that the most appropriate Navier- 
Stokes formulations for attention in the area of applica- 
tions involving realistic Reynolds numbers and reasonably 
complex geometries, are the time-dependent RANS and 
PNS formulations. The remainder of this survey will be 
solely devoted to these areas with primary emphasis on 
the time-dependent RANS formulation. For an indepth 
technical discussion of most of the CFD formulations just 
discussed, the interested reader is referred to Ref. 12. 


Previous Survey Work 

Several interesting surveys on the subject of Compu- 
tational Aerodynamics have been produced recently. Two 
particularly important articles are from Chapman 11 ’ 13 
where the entire field of Computational Aerodynamics is 
reviewed; pacing items including turbulence modeling, 
three-dimensional grid generation, algorithm development, 
and computer speed limitations are discussed in detail; and 
computational requirements and projected computer capa- 
bilities are examined. More recent surveys along the same 
line are by Kutler. 10,14 The same critical pacing items are 
examined and updated. 

Other pertinent reviews in the area of CFD are by 
Shang 15 where the RANS formulation is discussed for var- 
ious geometric applications, Mehta and Lomax 16 where 
the RANS formulation for transonic flow is discussed with 
special emphasis on algorithmic details and aspects of tur- 
bulent physics, Chaussee 17 where the PNS approach is 
described with special emphasis on applications, Orszag 
and Israeli 18 where early work associated with the solu- 
tion of the incompressible Navier-Stokes equations is re- 
viewed, and Turkel 19 where computational physics is re- 


viewed with emphasis on algorithms and convergence ac- 
celeration. Additional areas related to the numerical so- 
lution of the Navier-Stokes equations are grid generation, 
which is reviewed in detail by Thompson, 20 and turbu- 
lence modeling, which is described with particular empha- 
sis on three-dimensional flows involving curvature, separa- 
tion, and vorticies, by Lakshminarayana. 21 

An important aspect associated with CFD research is 
the area of experimental fluids research. A good review 
in this area is given by Marvin 22 where CFD code verifi- 
cation experiments are described and in Ref. 2 where the 
projected influence of CFD on experimental facilities is dis- 
cussed. Of course, the influence of the experimental world 
on CFD has been profound. The opposite is also true, CFD 
has begun to influence the experimental world as well. The 
CFD field has generated the need to conduct a whole new 
series of experiments to validate CFD code results. 22 These 
experiments generally take advantage of the latest develop- 
ments in sophisticated measurement techniques and gener- 
ally get more involved with fundamental flow phenomena. 
Having intimate CFD/experimental interaction is an im- 
portant feature and must remain as an essential ingredient 
in CFD research, especially in areas involving unexplained 
flow phenomena such as shock-induced or angle-of-attack 
induced flow separation. 

Before proceeding with the main emphasis of this re- 
view, an additional statement regarding experimental val- 
idation of numerical solutions is in order. A typical steady 
three-dimensional Navier-Stokes solution produced today 
consists of 100,000-200,000 grid points each with 5 depen- 
dent variables resulting in a total data base of perhaps one 
million words for each set of free-stream conditions. Cor- 
responding experimental data bases, not counting qualita- 
tive flow visualization results, generally consist of perhaps 
only a few hundred data points for each set of free-stream 
conditions. These data are usually in the form of surface 
pressures. This general lack of extensive data, including 
interior flow field data, makes thorough validation of CFD 
results difficult. If the experimental and computational 
surface pressures agree, this by no means requires other 
aspects of the flow field to agree. If the surface pressures 
disagree, the physical reasons for this disagreement are 
usually not apparent. Thus, it is imperative that sufficient 
detail associated with CFD validation experiments be ob- 
tained to fully evaluate numerical results and to provide 
the maximun understanding of basic fluid flow phenomena. 

Status of RANS/PNS CFD 

This survey is divided into several different areas of 
application as follows: incompressible flow, transonic flow, 
supersonic/hypersonic flow, and high-alpha/vortex-domi- 
nated flow. The f rst two areas will be limited to time- 
dependent RANS applications and the latter two will in- 
clude both time-dependent RANS and PNS applications. 
All the research surveyed will be restricted to three spa- 
tial dimensions and most of the applications will be asso- 
ciated with steady flow. First, the discussion centers on 
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the numerical solution of the Navier-Stokes equationsfor More recent work in the area of INS solvers has been 


incompressible flow. 


Incompressible Flow Applications 

Numerical solution of the incompressible Navier- 
Stokes (INS) equations in three dimensions has received 
much attention during the past several years. In the study 
of incompressible flow algorithms, it is interesting to study 
compressible algorithms and the connections between the 
two fields. The overwhelming majority of all steady com- 
pressible algorithms applied to flows that are partially 
or totally subsonic outside the boundary layer, utilize a 
time-accurate or time-like iteration to asymptotically ap- 
proach steady-state solutions with increasing time. The 
application of this approach to incompressible problems 
is difficult because of the basic difference in incompress- 
ible and compressible continuity equations. The incom- 
pressible continuity equation is the familiar condition that 
requires the divergence of the velocity vector to vanish. 
Thus, there is no time-dependent term as with the com- 
pressible continuity equation. Also associated with this 

problem is the lack of an explicit relation for the pres- 
sure. These two features cause the primary difference be- 
tween incompressible and compressible algorithms. Note: 
Certain two-dimensional algorithms, for example, stream- 
function/vorticity schemes, deal with these difficulties 
quite nicely but are difficult to extend to three dimensions, 
and therefore, will not be discussed further. 

Early numerical work in the incompressible flow area 
was presented by Harlow and Welch, 23 Chorin, 24 and 
Patankar and Spalding 25 where three basically different 
approaches were presented. In the first approach, a Pois- 
son equation for the pressure was derived by taking the 
divergence of the momentum equation. This equation is 
then solved iteratively using the velocity field at time level 
n to predict the pressure such that the continuity equation 
is satisfied at time level n+1. This process is time accu- 
rate and can be continued until a steady-state solution is 
reached. Needless to say, because of the additional itera- 
tion required to solve the pressure Poisson equation, this 
method is computationally expensive. 

In the second approach, the method of artificial com- 
pressibility was introduced. An artificial time-dependent 
pressure term is added to the steady continuity equation 
which allows the resulting equations to be solved using 
standard time-dependent compressible algorithms. As the 
solution converges, the time-dependent pressure term is 
driven to zero and the flow approaches the incompress- 
ible limit. In the last approach, a space marching proce- 
dure was developed by reducing the incompressible Navier- 
Stokes equations to a set of parabolic equations. The 
pressures are determined by first computing the velocities 
from estimated pressures, then a correction is applied to 
the pressure so as to satisfy the continuity equation. Of 
course, the latter approach does not work successfully if 
streamwise flow separation occurs. 


presented in Refs. 26-32 where three-dimensional ducts 
with curvature were considered for a wide variety of flow 
conditions, Himeno et al. 33 where the wake structure 
downstream of various bluff bodies was computed using 
the INS method described by Kawamura and Kuwahara, 34 
Kwak et al. 35 where a variety of applications were re- 
ported including a rectangular duct flow with a 45 deg 
bend, Gorski et al. 36 where turbulent corner flows were 
presented using two different turbulence models, and Refs. 
37-41 where various internal elements of the Space Shuttle 
Main Engine were analyzed. The algorithm used in Refs. 
37-41 was the method of pseduo compressibility devised 
by Chang and Kwak. 42 This method is similar to Chorin’s 
artificial compressibility method and is an outgrowth of 
the work of Steger and Kutler. 43 A recent upgrade to this 
algorithm was presented by Rogers et al. 44 where the orig- 
inal block-matrix ADI algorithm has been diagonalized to 
yield a scalar-matrix ADI algorithm which is about three 
times faster. 

A typical incompressible flow result showing particle 
traces downstream of a circular-cylinder post mounted on 
a flat plate is displayed in Fig. 5 (taken from Ref. 40). 
This calculation is laminar with a Reynolds number of 
1000. Note the formation of a saddle point of separation 
forward of the post and a horseshoe vortex downstream of 
the post. The secondary flow in front of the post wraps 
around the post and forms a “tornado-like” vortical flow 
in the wake of the post, which is strikingly different than 
the flow downstream of a simple post in the absence of 
the flat plate. An oil-flow photograph appearing in Ref. 
40 provides qualitative validation for this type of vortical 
wake flow. 



Fig. 5. Computed three-dimensional particle paths around 
and downstream of a circular-cylinder post mount- 
ed on a flat plate, incompressible flow, Re = 1000, 
Kwak et al. 40 
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Transonic Flow Applications 

The transonic flow area has perhaps received the wid- 
est attention during the past five years in the three-dimen- 
sional RANS research community. In this area there are 
a whole host of unsolved and interesting research topics 
that have practical applications for both commercial and 
military circles. A few of these include shock /bound ary- 
layer interation including boundary-layer separation ef- 
fects, wing/fuselage interference, performance prediction 
especially drag prediction, afterbody boundary-layer sep- 
aration, and helicopter rotor-tip transonic performance. 

Literally all transonic numerical methods for the 
Navier-Stokes equations utilize time- dependent formula- 
tions to obtain steady-state solutions asymptotically with 
increasing time. The time-integration algorithms used to 
implement this strategy vary from application to applica- 
tion. Most schemes utilize implicit methods, for example, 
the Beam-Warming, 45 Briley-McDonald, 46 or MacCor- 
mack 47 methods, or algorithms with suitable convergence 
acceleration, for example, multigrid. This theoretically 
allows large time steps to be taken, and therefore, pro- 
vides rapid approach to the desired asymptotic solution. 
This philosophy works well but requires several hundred 
to several thousand iterations to produce acceptable lev- 
els of “convergence.” The number of iterations generally 
increases as the mesh is refined and as the physics of the 
problem becomes more difficult, that is, stronger shock 
waves, larger separations, etc. Time-step limitations due 
to nonlinear numerical instabilities are common. This al- 
lows explicit time-integration methods, for example, 
MacCormack, 48 which generally have quite severe time- 
step limitations but which are cheaper per iteration and 
easier to code, to be appropriate alternatives for some ap- 
plications. 

The problem which has received the most attention in 
the transonic flow area is wing flow-field prediction. The 
first calculation of this type was presented by Mansour 49 
where the RANS formulation was solved for the flow field 
about “Wing C” (defined in Ref. 50). The numerical 
scheme used in this study was a fully-implicit, ADI-like 
scheme developed by Beam and Warming 45 and Briley and 
McDonald 46 and first demonstrated in three dimensions by 
Pulliam and Steger. 51 The grid used in this study, while 
having an efficient topology, was coarse and produced only 
fair agreement with experiment. Another similar calcula- 
tion for transonic wing geometries was presented by Agar- 
wal and Deese. 52 In this study the explicit Runge-Kutta 
scheme of Jameson et al. 53 was extended to handle the 
RANS equations for the transonic flow about the ONERA- 
M6 wing (defined in Ref. 54). The grid was again coarse, 
but for the simple calculations reported (attached and non- 
lifting), the experimental/computational agreement was 
good. 

More recent transonic RANS calculations for wing 
geometries have been given by Holst et al., 55 Flores, 56 
Kaynak et al., 57,58 and Srinivasan et al. 59 These studies 
have utilized a computer program called TNS (Transonic 
Navier-Stokes) which uses the diagonalized-implicit algo- 
rithm of Pulliam and Chaussee 60 coupled with a zonal grid 


arrangement to solve the thin-layer form of the RANS 
equations. In addition, the Baldwin-Lomax turbulence 
model 6 was used in each case. In the TNS approach for 
grid zones near the wing surface, the RANS equations are 
solved and for zones away from the wing surface, where 
viscous terms are not important, the Euler equations are 
solved. In Holst et al. 55 solutions are presented for a low- 
aspect-ratio wing derived from the NACA 0012 airfoil sec- 
tion. Comparisons with the experimental data of Lock- 
man and Seegmiller 61 for this wing at various transonic 
flow conditions including cases with shock/boundary-layer 
separation are presented and discussed. A significant as- 
pect of this set of calculations is that the wind tunnel walls 
are very important and must be modeled to obtain good 
agreement with experiment. Generally good agreement 
was obtained when the walls were modeled and the flow 
was attached. For separated cases the agreement was not 
good, and this was attributed to the lack of proper grid 
refinement and inappropriate turbulence modeling. 

Computed pressures produced by Kaynak et al. 57,58 
for the previously discussed Wing C geometry are com- 
pared with experimental results from two sources, Hin- 
son and Burdges 50 and Keener, 62,63 in Fig. 6. The free- 
stream Mach number for this set of results is 0.85, the 
angle of attack is 5 deg, and the Reynolds number is 8 
million. The computed and experimental results are in 
good agreement; as much scatter between the two exper- 
imental results exists as between the computed and ex- 
perimental results. Figure 7 shows the computed particle 
paths on the wing surface for the same case. Note the 
small shock-induced separation near the wing tip which is 
generally in good agreement with the corresponding ex- 
perimentally produced oil flow pattern. Detailed compu- 
tational/experimental comparisons for other cases are also 
presented in Kaynak et al. 57,58 

The computed results presented by Srinivasan et al. 59 
also deal with the Wing C geometry as well as various heli- 
copter rotor geometries. The main emphasis in this study 
is on tip vortex formation. With this purpose in mind spe- 
cial grid clustering at the wing tip was used. Computed 
particle paths on the surface of an ONERA wing designed 
to represent a typical helicopter rotor blade, are shown in 
Fig. 8. For this calculation the free-stream Mach number 
was 0.85, the angle of attack was 5 deg, and the Reynolds 
number was 8.5 million. A strong transonic shock wave 
near mid chord causes a large region of separated flow as 
shown in Fig. 8. An outward spiralling vortex emanates 
from the separated flow region near the tip. This can be 
seen in Figs. 9 and 10, which show standard and close-up 
views, respectively, of computed three-dimensional parti- 
cle paths about the tip. As explained by Srinivasan et 
al., 59 the cross flow created by the tip sweep enables the 
fluid particles released in the vincinity of the tip to first 
braid and then roll up and lift off of the surface. 

Other notable transonic wing calculations are pre- 
sented by Obayashi and Fujii 64 and Fujii and Obayashi. 65 
In the later case high-aspect-ratio transport-type-wing 
flow fields are obtained using a new LU-ADI algorithm. 
This algorithm utilizes the flux vector splitting techniau*' 
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Fig. 9. Computed three-dimensional particle paths around 

a swept-tip helicopter blade, Moo = 0.85, a = 5°, 
Re = 8.5xl0 6 , Srinivasan et al. 59 



Fig. 10. Blow-up view of Fig. 9 showing tip flow-field de- 
tail, Srinivasan et al. 59 


of Steger and Warming 66 and the diagonally dominate 
factorization of Lombard et al. 67 to decompose the usual 
block tridiagonal system into the product of lower and up- 
per scalar bidiagonal matricies. In this study flow fields 
with strong shock/boundary-layer interation and the re- 
sulting massive boundary layer separation are computed 
for a series of cases. Good correlations with experimen- 
tally measured pressures are presented. 

The LU-ADI transonic wing method just discussed 
has also been extended to handle wing/fuselage geometries 
by Fujii and Obayashi. 68 A typical result from this study is 
shown in Figs. 11 and 12. This case invloves a high- aspect- 
ratio wing low- mounted on a transport-type fuselage. The 
Mach number, angle of attack, and Reynolds number for 
this solution are 0.82, 4.0 deg, and 1.67 million, respec- 
tively. Computed and experimental pressure distributions 
are compared in Fig. 11 and computed surface particle 
paths are shown in Fig. 12. Note the shock-induced sepa- 






Fig. 11. Pressure coefficient comparisons for a high-aspect- 
ratio wing/fuselage configuration, Moo = 0.82, 
a = 4°, Re = 1.67xl0 6 , Fujii and Obayashi. 68 



Fig. 12. Computed surface particle paths for a high-aspect- 
ratio wing-fuselage configuration, Moo = 0.82, 
a = 4°, Re = 1.67xl0 6 , Fujii and Obayashi. 68 


ration at about mid chord and the wing-fuselage/juncture 
separation which is highlighted by a large spiral node on 
the aft-wing-root surface. This set of calculations is also 
quite unique due to the level of grid refinement used. For 
these calculations a grid of 119 x 71 x 92 = 777,308 points 
was used and represents the finest RANS grid with a pub- 
lished solution to date. 
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Additional RANS calculations have been published 
for a wide variety of other types of transonic applications. 
These studies include the simulation of nacelle/inlet flow 
fields by Vadyak, 69 calculations through swept-blade cas- 
cades by Swisshelm and Johnson, 70 the simulation of fore- 
body flow fields by Chaussee et al. 71 and Cosner, 72 the 
simulation of three-dimensional boattail afterbody flows 
by Deiwert, 73 and computations to determine the wake 
structure downstream of a three-dimensional turret by 
Purohit et al. 74,75 Another area with a lot of RANS re- 
search activity is transonic projectile aerodynamics. In 
this field the prediction of lift, drag, and magnus forces is 
very important and is intimately related to viscous effects. 
Interesting studies in this area include Nietubicz et al. 76 
where the Magnus forces associated with a spinning projec- 
tile were first computed, Nietubicz et al. 77 where detailed 
computations to predict various aerodynamic coefficients 
were presented, and Sahu et al. 78 where the projectile base 
flow field was computed at transonic speeds. This group 
of studies, ranging from inlets to projectiles, used a variety 
of different integration schemes ranging from the classical 
explicit MacCormack 48 scheme to all types of implicit and 
multigrid schemes. The turbulence models used are largely 
dominated by the Baldwin-Lomax model 6 with all of the 
models being in the algebraic category. The more com- 
plicated one- and two-equation models have not generally 
found their way into three-dimensional transonic Navier- 
Stokes applications. 


Supersonic/Hypersonic Flow Applications 

Supersonic flow field solutions are somewhat easier to 
obtain relative to transonic flow fields. This is because the 
physics of supersonic flow typically allows smaller compu- 
tational domains, the supersonic equation type is consis- 
tently hyperbolic (at least for inviscid regions) while for 
transonic flows the type is mixed between hyperbolic and 
elliptic, and the shock/boundary-layer interactions asso- 
ciated with supersonic flows are not as complicated as 
those associated with transonic flows. In addition, as al- 
ready discussed, the supersonic free-stream characteristic 
allows the introduction of marching algorithms, that is, 
the PNS formulation. Hypersonic flows, while enjoying 
the same advantages of supersonic flow, contain other com- 
plications that even today have not been fully evaluated. 
These include more complicated transition and turbulence 
physics, real gas effects, merged shock and boundary lay- 
ers, stronger shock waves, and radiation effects. 

The field of supersonic CFD utilizing both the time- 
dependent Navier- Stokes and PNS formulations for com- 
puting three-dimensional applications has been widely re- 
searched in the past 6-8 years. Early pioneering Navier- 
Stokes calculations in this area include Li 79 where lami- 
nar flow separation on blunt flared cones at angle of at- 
tack was studied, Holst and Tannehill 80 where viscous 
blunt body flows with an impinging shock were studied, 
Shang and Hankey 81 and Hung and MacCormack 82 where 
the flow in a three-dimensional compression corner was 
solved, Hung 83 where the flow over an inclined body of 
revolution was obtained, Hung 84 where the impingement 


of an oblique shock wave on a cylinder was investigated, 
Knight 85 and Hung and MacCormack 86 where an oblique- 
shock/boundary-layer interaction flow field was simulated, 
and Hung and Kordulla 87 and Hung and Buning 88 where 
the flow about a blunt-fin/flat-plate interaction was com- 
puted. The Mach number for these calculations ranged 
from low supersonic to hypersonic with the maximum be- 
ing 12.5. About half of the calculations, especially the 
eariler ones, were laminar, and the other half were tur- 
bulent with a variation of the Baldwin-Lomax turbulence 
model used in all cases. 

As seen in the previous examples, turbulence model- 
ing associated with the three-dimensional RANS formula- 
tion has not reached a very sophisticated level. For at- 
tached flow cases, existing algebraic, or “zero-equation” 
models as they are formally called, seem to do a good job 
for a wide variety of circumstances. When the flow expe- 
riences moderate or massive separation, these models do 
a poor job. For example, the most commonly used tur- 
bulence model is the Baldwin-Lomax model, which tends 
to underpredict the level and height of reversed flow in 
separation regions. This is probably due to an effective 
overprediction of dissipation in reqions of separated flow. 
More sophisticated models of the so-called one- or two- 
equation variety seem to improve the correlation in some 
cases while making it worse in others. Two difficulties as- 
sociated with their use are additional complication, that 
is, new differential equations have to be solved, and addi- 
tional numerical stiffness which causes the need for more 
flow solver iterations. Because of this lack of a clear di- 
rection and because of the additional complication, move- 
ment to use these more sophisticated turbulence models in 
three-dimensional applications has been slow. 

Additional supersonic Navier-Stokes calculations are 
given by Horstman et al. 89 for swept compression corners 
and axisymmetric bodies with either skewed or segmented 
conical flares; Horstman 90 for swept compression corner 
flows; and Knight, 91 Knight et al., 92 and Horstman 93 for 
swept-shock- wave/boundary- layer interaction studies. All 
of these three-dimensional computations involved Mach 3 
flows with significant shock-induced boundary layer sepa- 
ration. In addition, most of these studies utilize some form 
of the more sophisticated two-equation k-e eddy-viscosity 
turbulence model. In Ref. 92 computations with both the 
Baldwin-Lomax algebraic model and the Jones-Launder 94 
two-equation model were performed. The Jones-Launder 
computations also utilized Viegas-Rubesin 95 wall function 
boundary conditions which permit the use of coarse grids 
in the lower portion of the boundary layer. This reduces 
the amount of computer time required for a converged so- 
lution. Results from Ref. 92 indicate that, at least for 
the swept-shock/boundary- layer interaction flows studied, 
the overall structure of the three-dimensional flow was in- 
sensitive to the turbulence model, except within a small 
portion of the boundary layer near the surface. 

The difficulties associated with obtaining good turbu- 
lence models for complex applications exists for all speed 
regimes. However, the problems associated with the hy- 
personic flow regime are the most severe. At hypersonic 
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speeds certain incompressible flow assumptions used in for- 
mulating turbulence models, which are still valid up to 
Mach numbers of about 5, start to break down. In addi- 
tion, hypersonic flows, especially for Mach numbers well 
above 10, will have chemical reactions which may inter- 
act with the turbulence. For moderate to high altitudes, 
hypersonic flows will be transitional, further complicating 
the situation. It is clear that much work in the area of tur- 
bulence modeling in general, and for the hypersonic regime 
in particular, will be required in the near future. 

Flow field computations about more complicated ge- 
ometries have been explored in the supersonic/hypersonic 
flow regime. Examples of time-dependent RANS calcula- 
tions in this area include Shang 96 ’ 97 where the hypersonic 
wing-fuselage interference problem was studied with two 
different grid topologies, Kumar 98 where the supersonic 
flow upstream and through a two-strut scramjet inlet was 
computed, Howlett and Hunter 99 where the supersonic 
flow around and inside an inlet with spiked centerbody 
was simulated, Rizk and Ben-Shmuel 100 where the low- 
supersonic flow field about the space shuttle orbiter was 
investigated, and Li 101 where the hypersonic flow about an 
aerobrake body was considered. Another significant study 
in this area is the work of Shang and Scherr. 102 In this 
paper the entire flow field over the X-24C Lifting Body at 
a Mach number of 5.95, an angle of attack of 6 deg, and 
a Reynolds number of 16.4 million/m was obtained with 
a time-dependent RANS approach. The computational 
grid was extremely fine, consisting of 475,200 points, and 
represented the finest-grid RANS solution published up to 
1985. The explicit MacCormack numerical scheme was 
used and the turbulence was modeled with the Baldwin- 
Lomax model. A typical set of results from Ref. 102 is 
shown in Figs. 13-15. Figure 13 shows the surface stream- 
line pattern over the lifting body and demonstrates the 
level of complication that can be captured with the RANS 
approach when fine grids are used. Pressure distributions 
are compared with experiment in Figs. 14 and 15. Figure 
14 shows comparisons over the fore part of the lifting body 
and Fig. 15 comparisons over the aft part. For all stations 
the agreement is excellent. 

An additional complex-geometry study performed by 
Deiwert and Rothmund 103 involved the computation of 
flow about a three-dimensional body of revolution at angle 
of attack. This configuration consisted of a cone-cylinder 
forebody with a conical afterbody and contained a cen- 
tered, supersonic propulsive jet. Surface streamlines on 
the afterbody and density contours on the bilateral plane 
of symmetry are displayed in Fig. 16. For this compu- 
tation the free-stream Mach number was 2.0, the jet exit 
Mach number was 2.5, the angle of attack was 6 deg, and 
the jet-to-free-stream pressure ratio was 3.0. The flow was 
turbulent with a Reynolds number based on cylinder di- 
ameter of 1.5 million, and turbulence was modeled with 
the Baldwin-Lomax model. This particular set of calcula- 
tions is unique because of the complex data-base handling 
scheme utilized. The total data base associated with the 
216,000 point grid was about 5 million words and was di- 
vided into a series of blocks each consisting of 8000 grid 
points. Manipulations of these blocks were performed us- 



Fig. 13. Computed surface streamlines over the X-24C 
Lifting Body, Moo = 5.95, a = 6°, Re/m = 16.4 
xlO 6 , Shang and Scherr. 102 



NORMALIZED ARC LENGTH 


Fig. 14. Forebody pressure comparisons for the solution of 
Fig. 13, Shang and Scherr. 102 

ing the block data-base handling scheme described by Pul- 
liam and Lomax. 104 With this approach very large grids 
were accomodated on a computer with a limited main 
memory for a flow field of considerable complexity. 

The level of geometric complication, exhibited by the 
results from the preceeding paragraph and the section on 
transonic flow, has advanced rapidly during the past sev- 
eral years. However, a significant amount of work still 
remains. This area, along with turbulence modeling, was 
identified by Chapman 11 * 13 as one of the chief items pac- 
ing further development of computational aerodynamics. 
The key feature required to mature this area for three- 
dimensional applications is automation. Elements of this 
automation process must include surface grid generation, 
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SHANG ANDSCHERR 



Fig. 15. Afterbody pressure comparisons for the solution 
of Fig. 13, Shang and Scherr. 102 



Fig. 16. Afterbody flow detail: surface streamlines and 
density contours on the bilateral plane of symme- 
try, Moo = 2, Mj = 2.5, a = 6°, Pj/Poo = 3, 
Rep — 1.5xl0 6 , Deiwert and Rothmund. 103 


zoning algorithms, and interior grid generation with solu- 
tion adaptive capabilities. One possible alternative which 
may have a significant impact on the automation pro- 
cess is Artificial Intelligence as described by Andrews. 105 
Progress in individual applications is anticipated to pro- 
ceed at a fast rate, however, progress in automating the 
whole geometry and grid generation process will be slower 
to evolve. 

A large amount of work associated with viscous su- 
personic flows has been completed using the PNS formula- 
tion. This includes the pioneering work on cones and blunt 
cones at angle of attack by Lin and Rubin, 106 Lubbard 
and Helliwell, 107 and Agarwal and Rakich; 108 the work on 
jets by McDonald and Briley; 109 and the work on ogive- 
cylinders by Rakich et al. 110 Applications of the PNS for- 
mulation to supersonic projectile aerodynamics have been 
presented by Schiff and Sturek 111 and Sturek and Schiff 112 


using the formulation described in Schiff and Steger. 113 
Sensitive magnus force calculations are reported and de- 
scribed in the Ref. 112. 

Additional results utilizing the PNS formulation for 
more complicated geometries have been presented by Rai 
and Chaussee 114 and Rai et al. 115 for finned missies and 
projectiles, by Kaul and Chaussee 116 and Chaussee et 
al. 117 for the X-24C lifting body, by Chaussee et al. 118 for 
the Space Shuttle Orbiter, and by Chaussee et al. 119 and 
Wai et al. 120 for a generic supersonic-cruise fighter con- 
figuration. A typical result for the Space Shuttle Orbiter 
study of Ref. 118 is shown in Fig. 17. The free-stream 
Mach number for this calculation was 7.9, the angle of 
attack was 25 deg, and the Reynolds number was 2.4 mil- 
lion/m. Three-dimensional particle paths are displayed 
over the orbiter and show vortices forming on the leeside 
of the orbiter and at the wing/fuselage juncture. Another 
interesting PNS result is presented by Chaussee 17 for the 
blunt-nosed, biconic configuration displayed in Fig. 18. 
This configuration contained slices on both the windward 
and leeward sides and a control surface protrubing from 
the windward surface. The heat transfer distribution is 



Fig. 17. Computed three-dimensional particle paths about 
the Space Shuttle Orbiter configuration, PNS for- 
mulation, Mqo — 7.9, a — 25°, Re/m — 2.4xl0 6 , 
Chaussee. 17 




FLATCUT 



7° SLICE 
5° FLAP 


PLANFORM VIEW 


Fig. 18. Blunt biconic configuration with top and bottom 
slices and a flap control surface, Chaussee. 17 
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compared with experiment along two meridinal planes in 
Fig. 19. The Mach number for this calculation was 10, 
the angle of attack was 10 deg and the Reynolds num- 
ber was 0.3 million/m. Agreement is generally good with 
slight discrepancies caused by problems in predicting flow 
transition from laminar to turbulent. 


♦ EXPERIMENT, 0 = 90° PLANE 
4 EXPERIMENT. LEEWARD PLANE 



Fig. 19. Axial variation of the heat transfer over the biconic 
configuration in Fig. 18, PNS formulation, Moo = 
10, a = 10°, Re D = 8.3xl0 4 , Chaussee. 17 


High- Alpha/ Vortex-Dominated Flow Applications 

The field of high-alpha flows has recently gained a lot 
of momentum in the CFD community. Early work with 
PNS applications has been complemented with emerging 
activities in the time-dependent RANS area. High-alpha 
flows are more difficult to understand and to compute be- 
cause of several factors. First of all, associated with high- 
alpha flows are large vortex-dominated flow separations 
with complicated flow physics that are not well under- 
stood. Computing these flow fields accurately means pre- 
dicting vortex formation in the boundary layer, the lift-off 
process, the vorticity transport process, and the vortex- 
vortex interaction process. Because a significant amount 
of flow detail can exist out of the boundary layer, provid- 
ing the proper level of grid resolution in the proper places 
is more complicated. Secondly, high-alpha flows are gener- 
ally unsteady in nature and can be wildly asymmetric even 


when the boundary conditions of the flow are essentially 
symmetric. Turbulence modeling is another difficulty as- 
sociated with computing vortex-dominated flows. What 
turbulence models are appropriate for these flows, is a typ- 
ically asked question which is largely unanswered to date. 
A good survey for high-alpha, vortex-dominated flows is 
given by Newsome and Kandil. 121 

Much of the work on the PNS fromulation presented 
above in the section on supersonic/hypersonic flow is inti- 
mately related to the field of high- alpha flow field studies. 
For high-alpha flows which are supersonic and which con- 
tain no streamwise separations, utilization of the PNS ap- 
proach may be the most appropriate and certainly is more 
economical than the time-dependent RANS approach. 
Two examples that utilize this philosophy to study, in 
particular, high-alpha physics are given by Degani and 
Schiff 122 and Vigneron et al. 123 The results of Ref. 122 
are interesting in that the standard Baldwin-Lomax tur- 
bulence model was modified for high angle-of-attack flows 
involving vortex lift off. The boundary layer thickness is 
computed from the vorticity distribution in the standard 
Baldwin-Lomax model. This quantity is inappropriately 
determined when vorticies lift off the surface. The Degani- 
Schiff modification fixes this difficulty and produces good 
agreement with experiment even for vortex lift off. 

Computational work in the area of high-alpha flows 
has more recently centered on the use of the time-depen- 
dent RANS approach. Interesting examples in this area in- 
clude Pulliam and Steger, 51 Kordulla et al., 124 Obayashi 
et al., 125 Aki and Yamada, 126 and Ying et al. 127 where so- 
lutions over hemisphere-cylinder configurations have been 
studied in the transonic and low-supersonic regime; Fujii 
and Kutler, 128,129 Rizzetta and Shang, 130 Pan and 
Pulliam, 131 Mehta, 132 and Chaderjian 133 where high- 
alpha flows over wings, delta-wings, yawed wings, and 
wing-strake combinations have been studied; and New- 
some and Adams 134 where the vortical flow over elliptical 
bodies at high angles of attack have been studied. 

Results from Refs. 127 and 134 are particularly strik- 
ing and are displayed in Figs. 20 and 21, respectively. 
Figure 20 shows computed surface velocity vectors over a 



Fig. 20. Computed particle paths on the surface of a hemi-sphere-cylinder configuration, Moo 
4.45x10 s , Ying et al. 127 


1.2, a = 19°, Re D = 
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hemisphere-cylinder combination inserted into a Mach 1.2 
stream at an angle of attack of 19 deg and a Reynolds 
number based on the cylinder diameter of 445,000. This 
calculation was laminar and therefore no turbulence model 
was used. The surface velocity vector pattern, which is in 
good agreement with the experimental oil flow, was com- 
puted with a new upwind algorithm described in Ref. 127 


COMPUTATION 

□ EXPERIMENT 



Fig. 21. Pressure coefficient comparisons for an elliptical 
lifting body with angles of attack (a) and roll ( <t > ), 
Moo = 2.5, a = 20°, <t> = 45°, Re/m = 6.6xl0 6 , 
Newsome and Adams. 134 


on a fine grid consisting of (60 x 50 x 60 =) 180,000 points. 
Figure 21 shows computed pressure results from Ref. 134 
compared with experiment for the flow around an ellipti- 
cal body. The Mach number for this calculation was 2.5; 
the angles of attack and roll were 20 and 45 deg, respec- 
tively; and the Reynolds number was 6.6 million/m. The 
numerical grid for this calculation was quite fine consisting 
of (35 x 101 x 61 — ) 215,635 points. The explicit MacCor- 
mack method was used to integrate the RANS equations 
and a modified Baldwin-Lomax turbulence model was used 
to simulate turbulence. The excellent correlation of this 
calculation to experiment is obvious from Fig. 21. The 
reader is referred to additional interior flow field numer- 
ical/experimental comparisons, all in good agreement, in 
the original paper of Ref. 134. 


Another set of results dealing with the computation of 
vortical flows is the work of Reznick and Flores 135 and Flo- 
res et al. 136 In these studies the previously mentioned TNS 
computer code was upgraded to include more complicated 
geometries. In particular, a modified F-16A configuration 
with a fared inlet and the tail removed was simulated with 
a zonal grid topology consisting of 16 grid zones. Par- 


ticular emphasis in Ref. 135 was given to the study of 
the strake-leading-edge vortex. A comparison of computed 
wing pressures with experiment at a Mach number of 0.9, 
an angle of attack of 6 deg, and a Reynolds number of 6 
million is shown in Fig. 22. The two results are in generally 
good agreement. The discrepancies in shock location and 
leading-edge pressure levels are perhaps due to coarse grid 
effects or the fact that the experimental data contained a 
flow-through inlet that was not modeled computationally. 
Details of the vortex formation at the strake leading edge 
for a Mach number of 0.6 and an angle of attack of 10 deg 
are shown in Figs. 23 and 24. Figure 23 shows a view of 
the modified F-16A with selected three-dimensional parti- 
cle paths in the vicinity of the wing-root/strake juncture. 
Figure 24 shows a close-up view of the wing-root/strake 
juncture. From these two views, the formation of a vor- 
tex from both the strake and wing-root leading edges is 
obvious. The wing-root vortex draws the strake vortex 
beneath it and creates a strong spanwise flow toward the 
tip. This behavior is qualitatively displayed in experimen- 
tal flow visualizations described in Ref. 135. 


Concluding Remarks 


The numerical solution of the Navier-Stokes equations 
in three space dimensions has been reviewed. The areas of 
incompressible flow, transonic flow, supersonic/hypersonic 
flow, and high-alpha/ vortex-dominated flow have all been 
considered. Applications research in all of these areas is 
seen to be increasing at an exponential rate. In addition, 
the geometrical complication of applications is advancing 
but at a slower rate. Areas which need concentrated at- 
tention in the near future are listed as follows: 

1) Grid generation for complex geometries including 
surface grid generation, zoning algorithms, and solution 
adaptive grid concepts. 

2) New developments in turbulence models for ap- 
plications, especially those types of models that can be 
applied to obtain engineering answers for separated flow 
problems. 

3) More complete CFD validation experiments which 
will allow more complete validation of CFD results and 
a physical assessement of why disagreements occur when 
they do occur. 

4) More concentrated effort in the hypersonic flow re- 
gime including the development of sophisticated strong- 
shock capturing schemes, hypersonic turbulence models 
including the effects of transition, and the development of 
chemical equilibrium and nonequilibrium real-gas effects. 

5) Convergence acceleration techniques for time- 
asymptotic problems, especially new algorithms which will 
lend themselves to vectorization and be efficient for new 
more sophisticated turbulence models and the more com- 
plicated physics associated with hypersonic flow. 
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rj = 0.84 


Fig. 22. Pressure coefficient comparisons on the wing and fuselage of a modified F-16A configuration, = 0.9, a = 6.2 , 
Rec = 4.5xl0 6 , Reznick and Flores. 135 
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Fig. 23. Computed three-dimensional particle paths in the 
vicinity of a modified F-16A strake/wing juncture, 
Moo — 0-6, ex. = 10°, Reznick and Flores. 135 


Fig. 24. Blow-up view of the three-dimensional particle 
paths shown in Fig. 23, Reznick and Flores. 135 
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